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Abstract. We consider sequential random packing of cubes z + [0, 1]" with z G 
iZ" into the cube [0, 2]™ and the torus R n /2Z n as A — > oo. In the cube case [0, 2]" 
as N — ► oo the random cube packings thus obtained are reduced to a single cube 
with probability 1 — O (jf)- In the torus case the situation is different: for n < 2, 
sequential random cube packing yields cube tilings, but for n > 3 with strictly 
positive probability, one obtains non-extensible cube packings. 

So, we introduce the notion of combinatorial cube packing, which instead of de- 
pending on N depend on some parameters. We use use them to derive an expansion 
of the packing density in powers of . The explicit computation is done in the cube 
case. In the torus case, the situation is more complicate and we restrict ourselves 
to the case A — > oo of strictly positive probability. We prove the following results 
for torus combinatorial cube packings: 

• We give a general Cartesian product construction. 

• We prove that the number of parameters is at least "(" 2 +1 ) and we conjecture 
it to be at most 2" — 1. 

• We prove that cube packings with at least 2" — 3 cubes are extensible. 

• We find the minimal number of cubes in non-extensible cube packings for n 
odd and n < 6. 



1. Introduction 

Two cubes z + [0, l] n and z' + [0, l] n are non- overlapping if the relative interi- 
ors z+]0, l[ n and 2/+]0,l[ n are disjoints. A family of cubes (z % + [0, l] n )i<j< m with 
z l G jf^ n and N G Z>o is called a discrete cube packing if any two cubes are non- 
overlapping. We consider packing of cubes z + [0, l] n with z G into the cube 
[0, 2] n and the torus M. n /2Z n . In those two cases, two cubes z+ [0, l] n and z'+[0, l] n are 
non-overlapping if and only if there exist an index i G {1, . . . , n} such that Zi = z\ + 1 
(mod 2). A discrete cube packing is a tiling if the number of cubes is 2 n and it is 
non- extensible if it is maximal by inclusion with less than 2 n cubes. 

A sequential random cube packing consists of putting a cube z+ [0, 1]™ with z G -^T/ 1 
uniformly at random in the cube [0, 2] n or the torus K n /2Z, n until a maximal packing is 
obtained. Let us denote by M^(n), Mjj{n) the random variables of number of cubes of 
those non-extensible cube packings and by E(M^(n)), E(Mjj{n)) their expectation. 
We are interested in the limit — > oo and we prove that if N > 1 then 

(1) E{M%{n)) = (N k _i)k with U G T > and U M e Q 

fc=0 ^ ' 

In the cube case we prove that Ck(n) are polynomials of degree k, which we compute 



for k < 6 (see Theorem 3.3). In particular, Cq = 1, since as — > oo with probability 
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1 — O(jj), one cannot add any more cube after the first one. In the torus case the 
coefficients Tk(n) are no longer polynomials in the dimension n. The first coefficient 
T (n) = liniTv^oo E(Mj^(n)) is known only for n < 4 (see Table [l]). But we prove 
in Theorem 4.4 that if n > 3 then T (n) < 2™. This upper bound is related to the 
existence in dimension n > 3 of non-extensible torus cube packings (see Figure [2j 
Table [TJ Theorem 4.4 and Section [5]). 

Those results are derived using the notion of combinatorial cube packings which is 
introduced in Section |2j A combinatorial cube packing does not depend on N but 
instead on some parameters t{] to a cube or torus discrete cube packing CV, one can 
associate a combinatorial cube packing CV' = <f)(CV). Given a combinatorial cube 
packing CV the probability p(CV, N) of obtaining a discrete cube packing CV' with 
<p(CV') = CV is a fractional function of N. We say that CV is obtained with strictly 
positive probability if the limit hin^^oo p(CV, N) is strictly positive. 

In Section [3] the method of combinatorial cube packings is applied to the cube case 
and the polynomials Ck are computed for k < 6. In the torus case, the situation is 
more complicated and we restrict ourselves to the case of strictly positive probability, 
i.e. the limit case N — > oo. In Section [4] we consider a Cartesian product construction 
for continuous cube packings obtained with strictly positive probability. The related 
lamination construction is used to derive an upper bound on E(M^(n)) in Theorem 

m 

In Section [5j we consider properties of non-extensible combinatorial torus cube 
packings. Firstly, we prove in Theorem 5.1 that combinatorial cube packings with 
at least 2 n — 3 cubes are extensible to tilings. In Propositions 5J3 and 5.5[ we prove 
that non-extensible combinatorial torus cube packings obtained with strictly positive 
probability have at least "("+-*-) parameters and that this number is attained by a 
combinatorial cube packing with n + 1 cubes if n is odd. We conjecture that the 
number of parameters is at most 2" — 1 (see Conjecture 5.4). In Proposition 5.7 we 



prove that in dimension 6 the minimal number of cubes in non-extensible combinato- 
rial cube packings is 8 and that none of those cube packings is attained with strictly 



positive probability. In Proposition 5J3 we show that in dimension 3, 5, 7 and 9, there 
exist combinatorial cube tilings obtained with strictly positive probability and ra ^ n 2 +1 - > 
parameters. 

We now explain the origin of the model considered here. Palasti |Pa60j considered 
maximal packings obtained from random packings of cubes [0, l] n into [0,x] n . She 
conjectured that the expectation E(M x {n)) of the packing density M x {n) satisfies the 
limit 

(2) lim E ^)l = ft, . 

x— >oo x 

with f3 n — [3™. The value of f3\ is known since the work of Renyi [Re58j and in Penrose 
|Pe01j the limit ^ is proved to exist. Note that based on simulations it is expected 
that j3 n > ft™ and an experimental formula from simulations 

is known |BlSo82j . 

The Itoh Ueda model |ItUe83] is a variant of the above: one considers packing 
of cubes z + [0,2] n with z e Z n into [0,4] n . It is proved in |DIP051 [PoT)5l lPoT)3] 
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that the average number of cubes satisfies the inequality E(M% (n)) > (|) n and some 
computer estimate of the average density ±E{M$(n)) were obtained in |ItSo86j . In 
|DIP06j . we considered the torus case, similar questions to the one of this paper and 
a measure of regularity called second moment, which has no equivalent here. 

2. Combinatorial cube packings 

If z + [0, l] n C [0, 2} n and z = (z u . . . , z n ) G ±Z n then z { G {0, . . . , 1}. Take 
a discrete cube packing CV = (z l + [0, l] n )i<i<m of [0,2] n . For a given coordinate 
1 < j < n we set (j>(zj) = Uj with tij a parameter if < Zj < 1 and <f>(zj) = z* if 
Zj = or 1. If z % = (z\,zl,..., z l n ) then we set 0(V) = (<j>(z\), . . . , 4>{z l n )) and to CV 
we associate the combinatorial cube packing <fi(CV) = (4>(z l ) + [0, l] n )i<j< m . 

Take a torus discrete cube packing CV = (z % + [0, l] n )i<j< m with z % G ^Z n . For a 
given coordinate 1 < j < n we set = tfcj if 2] = ^ (mod 2) and = tk,j+l if 
2* = jj + 1 (mod 2) with tkj a parameter. Similarly, we set <t>(z % ) = (4>(z{), . . . , 4>(z l n )) 
and we define 4>{CV) = (4>(z l ) + [0, l] n )i<j< m the associated torus combinatorial cube 
packing. 

In the remainder of this paper we do not use the above parameters but instead 
renumber them into t±, . . . ,tjf. Without loss of generality, we will always assume that 
different coordinates have different parameters. Of course we can define combinatorial 
cube packing without using to discrete cube packings. In the cube case, the relevant 
cubes are of the form z + [0, 1]™ with zi = 0, 1 or some parameter t. In the torus case, 
the relevant cubes are of the form z + [0, l] n with Zi — t or t + 1 and t a parameter. 
Two cubes z* + [0, l] n and z* + [0, l] n are non-overlapping if there exist a coordinate 
j such that Zj = Zj +1 (mod 2). In the cube case this means that Zj = or 1 and 
Zj — 1 — Zj. In the torus case this means that Zj depends on the same parameter, 
say t, Zj, Zj = t or t + 1 and Zj ^ Zj. A combinatorial cube packing is then a family 
of such cubes with any two of them being non-overlapping. Notions of tilings and 
extensibility are defined as well. Moreover, a discrete cube packing is extensible if 
and only if its associated combinatorial cube packing is extensible. Denote by m(CV) 
the number of cubes of a combinatorial cube packing CV and by N(CV) its number of 
parameters. Denote by Comb c (n), Comb T (n), the set of combinatorial cube packings 
of [0, 2] n , respectively M n /2Z n . 

Given two combinatorial cube packings CV and CV' (either on cube or torus), 
we say that CV' is a subtype of CV if after assigning the parameter of CV to 0, 
1, or some parameter of CV' , we get CV' . So, necessarily m{CV') = m{CV) and 
N(CV') < N(CV) but the reverse implication is not true in general. A combinatorial 
cube packing is said to be maximal if it is not the subtype of any other combinatorial 
cube packing. Necessarily, a combinatorial cube packing CV is a subtype of at least 
one maximal combinatorial cube packing CV'. 

Given a combinatorial cube packing CV the number of discrete cube packings 
CV' such that <fi(CV') = CV is denoted by Nb(CV,N). In the cube case we have 
Nb(CV, N) = (N — \) N ( CV \ The torus case is more complex but it is still possible to 
write explicit formulas: denote by Nj(CV) the number of parameters which occurs in 
the j-th coordinate of CV. We then get: 

(3) Nb(CV, N) = U] =l U^if v) {2N - 2(k - 1)). 
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The asymptotic order of Nb(CP,N) is (2N) N{ - CV \ which shows that Nb(CP,N) > 
for N large enough. More specifically, Nj(CV) < 2 n so Nb(CP,N) > if iV > 2 n . 
Note that it is possible to have CP' a subtype of CP and Nb(CV', N) > Nb(CP, N) 
for small enough N. 

Denote by fjj(n) the minimal number of cubes of non-extensible discrete torus cube 
packings (z* + [0, l] n )i<j< m with z l G \ Denote by f^(n) the minimal number of 
cubes of non-extensible combinatorial torus cube packings. 

Proposition 2.1. For n > 1 we have hin^^oo fjf{n) = f^{n). 

Proof. A discrete cube packing CP is extensible if and only if <p(CP) is extensible. 
Thus f^(n) > fjj(n). Take CP a non-extensible combinatorial torus cube packing 
with the minimal number of cubes. By Formula ^ there exist No such that for 
N > N we have Nb(CP, N) > 0. The discrete cube packings CP' with <j>{CV) = CP 
are non-extensible. So, we have liniTv^oc fjf( n ) — f^( n )- D 
In the cube case we have for N >2 the equality f^(n) = 1. 

Two combinatorial cube packings CP and CP' are said to be equivalent if after 
a renumbering of the coordinates, parameters and cubes of CP one gets CP' . The 
automorphism group of a combinatorial cube packing is the group of equivalences of 
CP preserving it. Testing equivalences and computing stabilizers can be done using 
the program nauty [MKa05j . which is a graph theory program for testing whether 
two graphs are isomorphic or not and computing the automorphism group. The 
method is to associate to a given combinatorial cube packing CP a graph Gr{CP), 
which characterize isomorphism and automorphisms. The method used to find such 
a graph Gr(CP) are explained in the user manual of nauty and the corresponding 
programs are available from |Du07] . 

We now explain the sequential random cube packing. Given a discrete cube packing 
CP = (V + [0, l] n )i<i< m denote by Poss(CP) the set of cubes z + [0, l] n with z e ±Z n 
which do not overlap with CP. Every possible cube z + [0, l] n is selected with equal 
probability \p os h C p-\\ ■ The sequential random cube packing process is thus a process 
that add cubes until the discrete cube packing is non-extensible or is a tiling. 

Fix a combinatorial cube packing CP, N > 2 n and a discrete cube packing CP' 
such that 4>{CP') = CP. To any cube w + [0, l] n G Poss(CP') we associate the 
combinatorial cube packing CP W = 4>{CP' U {w + [0, l] n }). The set Poss{CP') is 
partitioned into classes Cli, . . . , Cl r with two cubes w + [0, l] n and w' + [0, l] n in 
the same class if CP W = CP W >. The combinatorial cube packing associated to Cli is 
denoted by CP{. The set {CPi, . . . ,CP r } of classes depends only on CP. If we had 
chosen some N < 2 n , then some of the preceding CPi might not have occurred. So, 
we have 

\Cl (N)\- Nb{CV - N) 
|CMiV)| - Nb(CP,N) 

and we can define the probability p(CP, CPi, N) of obtaining a discrete cube packing 
of combinatorial type CPi from a discrete cube packing of combinatorial type CP: 

p(C p CP N) = \<MW\ = NbjCP^N) 

n ' l) ) \Ch(N)\ + ■ ■ ■ + \Cl r (N)\ Nb{CVi,N) + --- + Nb(CV r ,N)' 

Given a combinatorial cube packing CP with m cubes a pathp = {CP , CP 1 , . . . , CP m } 
is a way of obtaining CP by adding one cube at a time starting from CP = and 
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Figure 1. Two 2-dimensional torus combinatorial cube tilings 



ending at CV m = CV. The probability to obtain CV along a path p is 

p(CV,p,N) = p(CV°,CV 1 ,N) x ■■■xp{CV m -\CV m ) N). 

The probability p(CV, N) to obtain CV is the sum over all the paths p leading to CV 
of p(CV,p, N). The probabilities p(CV,p, N) and p(CV,N) are fractional functions 
of N, which implies that the limit value p{CV, oo), p{CV,p, oo) and p(CV, CV' , oo) 
are well defined. 

As A goes to oo we have the asymptotic behavior 

|C7;(A)| ~ (2A) nfe * 

with n&j = N{CVi) ~ N{CV) the number of new parameters in CPj as compared with 
CV. Clearly as N goes to oo only the classes with the largest nbi have p(CV, CVi, oo) > 
0. If Ck is such a class then we get 

p(CV, CVi, oo) = ^ 

with r' the number of classes CU having the largest nbi and otherwise p(CV, CVi, oo) = 
0. Analogously, for a path p leading to CV we can define p{CV,p, oo) and p{CV , oo). 
We say that a combinatorial cube packing CV is obtained with strictly positive prob- 
ability if p(CV, oo) > that is for at least one path p we have p{CV,p, oo) > 0. For a 
path p = {CV°, CV 1 ,..., CV m } we have p{CV,p, oo) > if and only if every CV 1 has 
N{CV l ) maximal among all possible extensions from CV 1 ' 1 . This implies that each 
CV 1 is maximal, i.e. is not the subtype of another type. As a consequence, we can 
define a sequential random cube packing process for combinatorial cube packing CV 
obtained with strictly positive probability and compute their probability p(CV, oo). 

A combinatorial cube packing CV is said to have order k = ord(C'P) if p{CV , N) = 
(Af-i)fc f(N) with limAr^oo f(N) G M.* + . A combinatorial cube packing is of order if 
and only if it is obtained with strictly positive probability. 

Let us denote by M^(n), Mjj(n) the random variables of number of cubes of those 
non-extensible cube packings and by E(M^(n)), E(Mjf{n)) their expectation. From 
the preceding discussion we have 

E{M%{n)) = P( CV > N)m{CV) with U G {C, T}. 

CVeComb u (n) 

Denote by /J 0oo (n) the minimal number of cubes of non-extensible combinatorial 
torus cube packings obtained with strictly positive probability. 
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Figure 2. The 3-dimensional combinatorial cube packings obtained 
with strictly positive probability; two laminations over 2-dimensional 
cube tilings, the rod tiling and the smallest non-extensible cube packing 
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Table 1. Number of packings and tilings for the case N = oo and 
N = 2 (see |DIP06j ) 



In dimension 2 (see Figure [I]), there are three combinatorial cube tilings. One of 
them is attained with probability 0; it is a subtype of the remaining two which are 
equivalent and attained with probability |. By applying the random cube packing 
process and doing reduction by isomorphism, one obtains the 3-dimensional combi- 
natorial cube packings obtained with strictly positive probability (see Figure |2|). The 
non-extensible cube packing shown on this figure already occurs in |La00l |DIP06| . In 
dimension 4, the same enumeration method works (see Table [T]) but dimension 5 is 
computationally too difficult to enumerate. 



3. Discrete Random cube packings of the cube 

We compute here the polynomials Cfc(n) occurring in Equation ([T]) for k < 6. We 
compute the first three polynomials by an elementary method. 

Lemma 3.1. Put the cube z 1 + [0, l] n in [0,2] ra and write 

I — {i : z\ = or N} 

then do sequential random discrete cube packing. 

(i) The minimal number of cubes in the packing is |/| + 1. 
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(ii) The expected number of cubes in the packing is \I\ + 1 + O ■ 

Proof. Let us prove (i). If |/| = 0, then clearly one cannot insert any more cubes. 
We will assume |/| > 1 and do a reasoning by induction on |/|. If one puts another 
cube z 2 + [0, l] n , there should exist an index i E I such that \z 2 — z}\ = 1. Take an 
index j ^ i such that z 2 is or 1. The set of possibilities to add a subsequent cube 

-}. So, one can assume that for j ^ i, 



N- 



is larger if Zj G {0, 1} than if Zj G {^, . . . . , 
one has < z 2 < 1. This means that any cube z + [0, l] n in subsequent insertion 
should satisfy \z^ — z 2 \ = 1, i.e. z,i = z\. So, the sequential random cube packing 
can be done in one dimension less, starting with z' 1 — (zi, . . . , Zj_i, Zi+i, . . . , z n ). The 
induction hypothesis applies. Assertion (ii) follows easily by looking at the above 
process. For a given i the choice of z 2 with < z 2 < 1 for j ^ i is the one with 
probability 1 — (]v+i)- So, all neglected possibilities have probability O (jy^j) and 
with probability 1 — (jyVj) the number of cubes is the minimal possible. □ 
See below the 2-dimensional possibilities: 



5 



The random variable M^-(n) is the number of cubes in the obtained non-extensible 
cube-packing. E(M^(n)) is the expected number of cubes and E(M^(n) \ k) the 
expected number of cubes obtained by imposing the condition that the first cube 
z 1 + [0, l] n has \{i : z\ = or 1}| = k. 

Theorem 3.2. For any n > 1, we have 



E(M%(n)) = 1 + 



2n 4n(n — 1) 



+ 



N + l (N+l) 2 



+ 



N + l 



as N — > oo. 



Proof. If one chooses a vector z in {0, . . . ,N} n the probability that \{i 



or iV}| = k is {j^iY {^+1}" * (l)- Conditioning over k e {0,1,..., n}, one 
obtains 



(4) 



n 



k=0 1 



' N- l \ n - k (n 
N+l J U 



E{M%( n ) | fc). 



So one gets 



E(M°(n)) = (§+}) n E(M^n)\0)+n 1 ^- 1 (^r 1 



E(M%(n) | 1) 



+ n(n-l)^(gr«n)|2)+0(^) ! 
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Clearly E(M^{n) | 0) = 1 and E(Mg(n) | 1) = 1 + E(M^(n - 1)). By Lemma |[ 
E{M°{n) | 2) = 3 + O (^) . Then, one has E{M°{n)) = 1 + (^) and 

E{M°{n)) = (1-^)^^(1-^)^2 + 0^)) 



^ I (A/+1) 2 



V-J. ^ ((W+1) 2 ) } + Af+1 + ^ ((Af+l) 



1 + 



2/i 



+ °((AT+1 



AT+1 1 ^ ^(Af+1) 2 

Inserting this expression into E(M^(n)) and Formula Q one gets the result. □ 
So, we get Coin) = 1, Ci(n) = 2n and C 2 (n) = An{n — 2). In order to compute 
Cfe(n) in general we use methods similar to the ones of Section [2j Given a cube 
z + [0, l] n with 2^ G {0, jf, . . . , 1} we define a face of the cube [0, l] n in the following 
way: if Zj = or 1 then we set if>(zi) = or 1 whereas if < Zi < 1 we set if)(zi) = U 
with ti a parameter. When the parameters ti of the vector (ip(zi), . . . , ip(z n )) vary in 
]0, 1[ this vector describes a face of the cube [0, l] n , which we denote by i]){z). This 
construction was presented for the first time in |Po03l IPo05] . 

If F and F' are two faces of [0, l] n , then we say that F is a sub-face of F' and write 
F C F' if F is included in the closure of F' . A subcomplex of the hypercube [0, l] n is 
a set of faces, which contains all its sub-faces. If CP is a cube packing in [0, 2] n , then 
the vectors z such that z + [0, l] n is a cube which we can add to it are indexed by 
the faces of a subcomplex [0, l] n with the dimension giving the exponent of (N —l) k . 
The dimension of a complex is the highest dimension of its faces. Given a discrete 
cube packing CV, we have seen in Section [2] that the size of Poss(CV) depends only 
on the combinatorial type <p(CV). In the cube case which we consider in this section 
Poss(CV) itself depends only on the combinatorial type. 

Theorem 3.3. There exist polynomials Ck(n) of n with deg Ck = k such that for 
any n and N > 1 one has: 

C k {n) 



(N — l) k ' 



E{M c N {n)) = 

k=0 

The polynomials Ck(n) are given in Table^ 

Proof. The image ip(Poss(CV)) is an union of faces of [0, l] n , i.e. a subcomplex of 
the complex [0, l] n . Denote by dim(F) the dimension of a face F of the cube [0, l] n . 
Denote by Poss(F) the set of vectors z G {0, j^, . . . , l} n with ip(z) = F. we have the 
formula: 

\Poss(F)\ = (N - l) dim ( F ) and \Poss(CV)\ = ^(N - l) dim ( F ). 

F 

The cubes, whose corresponding face in [0, 1]™ have dimension dim(ijj(Poss(CV))) 
have the highest probability of being obtained. If one seeks the expansion of E(M^(n)) 
up to order k and if CV is of order ord(CP) then we need to compute the faces of 
ip(Poss{CV)) of dimension at least &\m(ip (Poss(CV))) — (k — ord(CP)). The proba- 
bilities are then obtained in the following way: 

(at — ly™^) 

( 5 ) p^ n ) = y rr " " ( N -i)«n*Fr 

^F'Gip(Poss(CT)) With dim(F')>dim(i>(Poss(CP)))-(k-ovd(CV))^ > 
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The enumeration algorithm is then the following: 



Input: Exponent k. 

Output: List C of all inequivalent combinatorial types of non-extensible cube 
packings CP with order at most k and their probabilities p(CV, N) with an error of 
p( i 

r^{0}. 

£ ^0 

while there is a CP £ T do 

t^t\ {cv} 

ip(Poss(CV)) <— the complex of all possibilities of adding a cube to CV 
T <— the faces of ip(Poss(CV)) of dimension at least 

d\m(i)(Poss{CV))) - (k - ord(CP)) 
if T = then 

if CV is equivalent to a CV' in C then 

j9(CP', iV) «- p(CP', AT) + p(CP, iV) 
else 

£ «- £ U {CP} 
end if 
else 

for C e T do 

CV new ^CVU{C} 
p(CV neW} N) ^ p(CV,N)p(C,N) 
if CVnew is equivalent to a CP' in T then 

p(CP', AT) <- p(CP', AT) + ^(CP^, AT) 
else 

r <- r u {cv new } 

end if 
end for 
end if 
end while 



Let us prove that the coefficients Ck(n) are polynomials in the dimension n. If C is 
the cube [0, l] n then the number of faces of codimension I is 2M?J, i.e. a polynomial 
in n of degree I. Suppose that a cube packing CP = (V + [0, l] n )i<j< m has < 
Zj < 1 for n' < j < n. Then all faces F of i/j(Poss(CV)) of maximal dimension 
d = dim(ip(Poss(CV))) have < Zj < 1 for n' < j < n and z & F. When one chooses 
a subface of F of dimension d — I, we have to choose some coordinates j to be equal 
to or 1. Denote by I' the number of such coordinates j with n' < j < n. There 
are 2 1 ' ( n+ ; r" ) choices and they are all equivalent. There are still I — I' choices to be 
made for j < n' — 1 but this number is finite so in all cases the faces of ip(Poss(CV)) 
of dimension at least d — I can be grouped in a finite number of classes with the 
size of the classes depending on n polynomially. Moreover, the number of classes of 
dimension d is finite so the term of higher order in the denominator of Equation ^ 
is constant and the coefficients of the expansion of p(F, N) are polynomial in n. □ 
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k 


\Comb% 




C k (n) 




1 

2 
3 

A 

-± 

5 
6 


1 
2 
3 

7 

18 
86 

1980 


1 

In 
4n(n - 2) 
|{28n 3 - 153n 2 + 149n} 

1 r °nifirj 4 ^lAlftn 3 1 W7IY1 -n 2 4n7°1n1 

7 {208724n 5 ' 3516724n 4 + 18627854n 3 35643809™ 2 + 20444915™} 
283X15*^113 13 17 {1929868729224214329703n 6 46928283796201160537385n 5 
+397379056595496330171955n 4 - 1442659974291080413770375n 3 
+2205275555952621337847422™ 2 - 1115911322466787143241320n} 



Table 2. The polynomials Ck(n). \Comb k \ is the number of types of 
combinatorial cube packings CP with ord(CP) < k 



4. Combinatorial torus cube packings and lamination construction 

Lemma 4.1. Let CP be a non- extensible combinatorial torus cube packing. 

(i) Every parameter t of CP occurs, which occurs as t also occurs as t + 1. 

(ii) Let Ci, . . . , Cfc be cubes of CP and C a cube which does not overlap with CP' = 
CP — {Ci, . . . , Cfc}. The number of parameters of C , which does not occur in CP' is 
at most k — 1. 

Proof, (i) Suppose that a parameter t of CP occurs as t but not as t + 1 in the 
coordinates of the cubes. Let C — z + [0, l] n be a cube having t in its j-th coordinate. 
If C = z'+[0, l] n is a cube of CP, then there exist a coordinate f such that 2y = Zf+1 
(mod 2). Necessarily f ^ j since t + 1 does not occur, so C + ej does not overlap 
with C as well and obviously C + ej does not overlap with C . 

(ii) Let C = z + [0, l] n be a cube which does not overlap with the cubes of CP' . 
Suppose that z has k coordinates i\ < ■ ■ ■ < ik such that their parameters ti,...,tk do 
not occur in CP' . If Cj = z^ + [0, l] n , then we fix z i] = z\. + 1 (mod 2) for 1 < j < k 
so that C does not overlap with CP. This contradicts the fact that CP is extensible 
so z has at most k — 1 parameters, which do not occur in CP' . □ 

Take two combinatorial torus cube packings CP = (z l + [0, l] n )i<i<m and CP' = 
{z n + [0, 1]"') i<j<m'- Denote by (z /lJ + [0, l] n )i<j< m ' with 1 < i < m m independent 
copies of CP'; that is every parameter t' k of z' 3 is replaced by a parameter t\ k in z' 1 ' 3 . 
One defines the combinatorial torus cube packing CP k CP' by 

{z\ z' iJ ) + [0, l} n+n ' for 1 < i < m and 1 < j < m. 

Denote by CP\ the 1-dimensional combinatorial packing formed by (t + [0, 1], t + 1 + 
[0, 1]). The combinatorial cube packings CP\ x CP\ and CP\ x (CP\ x CPi) are the 
ones on the left of Figure [T] and [2j respectively. Note that in general CP tx CP' is not 
isomorphic to CP' x CP. 

Theorem 4.2. Let CP and CP' be two combinatorial torus cube packings of dimension 
n and n' , respectively. 

(i) m(CP x CP') = m(CP)m(CP') and N(CP x CP') = N(CP) +m(CP)N(CP'). 

(ii) CP x CP' is extensible if and only if CP and CP' are extensible. 
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(iii) If CP and CP' are obtained with strictly positive probability and CP is non- 
extensible then CP x CP' is attained with strictly positive probability. 

(iv) One has /£(n + m) < fZ(n)fL( m ) and f> 0>oo (n + m) < /> ,oo(™)/>o,oc( m )- 

Proof. Denote by (z i + [0, l] n )i<i< m and by (z' j + [0, l] n ')i<j< m > the cubes of CP and 
CP' obtained in this order, i.e. first z 1 + [0, 1]™, then z 2 + [0, l] n and so on. Assertion 
(i) follows by simple counting. 

If CP, respectively CP' is extensible to CP U {C}, CP' U {C} then CP x CP' is 
extensible to (CP U {C7}) x CP', respectively CP x (CP' U {C}) and so extensible. 
Suppose now that CP and CP' are non-extensible and take a cube z+ [0, l] n+n with z 
expressed in terms of the parameters of CP x CP' . Then the cube (z\, . . . , z n ) + [0, l] n 
overlaps with one cube of CP, say z % + [0, l] n . Also (z n+ i, . . . , z n+n i) + [0, 1]" overlaps 
with one cube of CP' , say z' 3 + [0, l] n . So, z + [0, l] n+n overlaps with the cube 
(z\z' i,j ) + [0, l] n+n ' and CP x CP' is non-extensible, establishing (ii). 

A priori there is no simple relation between p(CPt< CP' , oo) and p(CP, oo), p(CP' , oo) 
But we will prove that if p(CP, oo) > 0, p(CP' , oo) > and CP is not extensible then 
p(CP x CP' , oo) > 0. That is, to prove (iii) we have to provide one path, among pos- 
sible many, in the random sequential cube packing process to obtain CP x CP' with 
strictly positive probability from some corresponding paths of CP and CP'. We first 
prove that we can obtain the cubes ((z\ z' hl ) + [0, l] n+n ')i<j< m with strictly positive 
probability in this order. Suppose that we add a cube z + [0, l] n+n after the cubes 
(z 1 ' , z' % ' ) + [0, l] n+n ' with i' < i. If we choose a coordinate k e {n + 1, . . . ,n + n'} 
such that Zfc = (z 1 ' , z n ' 1 )/ c + 1 for some i' < i then we still have to choose a coordinate 
for all other cubes. This is because all parameters in (z' h )i<i< m are distinct. So, 
we do not gain anything in terms of dimension by choosing k G {n + 1, . . . ,n + n'} 
and the choice (z % , z n ' ) has the same or higher dimension. So, we can get the cubes 
((z l , z' 1 ' ) + [0, l] n+n ')i<i< m with strictly positive probability. 

Suppose that we have the cubes (z l , z' lJ ) + [0, l] n+n ' for 1 < i < m and 1 < j < m' . 
We will prove by induction that we can add the cubes ((z\ z' l,m ° +1 ) + [0, l] n+n ')i<j< m . 
Denote by n' , < n' — 1 the dimension of choices in the combinatorial torus cube 

packing (z'° + [0, l] n ')i<i<m^ 

Let z + [0, l] n+n be a cube, which we want to add to the existing cube packing. 
Denote by S z the set of i such that z + [0, V\ n+n does not overlap with (z l , z n,J ) + 
[0, l] n+n on a coordinate k < n. The fact that z+ [0, l] n+n does not overlap with the 
cubes (z\ z nj ) + [0, l] n+n fixes n' —n'i coordinates of z. If i ^ i' then the parameters 

in z' 1 ' 3 and z n ,3 are different; this means that (n' — n' m ,)\S z \ components of z are 
determined. Therefore, since CP is non-extensible, we can use Lemma 4.1 (ii) and so 
get the following estimate on the dimension D of choices: 

D < W - [ri - n> m , Q )\S z \} + {\S z \-l} 
(6) < n' m ,-(n>-n> m , - 1){|^| - 1} 

< n' ,. 

We conclude that we cannot do better in terms of dimension than adding the cubes 
((z l , z' l,m ° +1 ) + [0, l} n+n ')i<i<m, which we do. So we have a path p with p(CP x 
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CV' ,p, oo) > which proves that CV x CV' is obtained with strictly positive proba- 
bility. 

(iv) follows immediately from (hi) and (ii). □ 
There exist cube packings CP, CV' obtained with strictly positive probability such 

that p(CV x CV' , oo) > 0, which shows that the hypothesis CV non-extensible is 

necessary in (iii). 

The third 3-dimensional cube packings of Figure |2j named rod packing has the 
cubes (h l + [0, l] 3 )i<i<s with the following h l : 

h 1 = ( h, t 2 , t 3 ) h 5 = ( t 6 + l, t 2 + l, t 5 + l) 

h 2 = ( ti + 1, U, h) h 6 = ( t x , t 2 , t 3 + l) 

h 3 = ( t 6 , h + 1, t 5 + 1) h 7 = ( t x + 1, t 2 , h + 1) 

h 4 = ( ti + 1, t 4 + l, *s) ^ 8 = ( *i, *2 + l, f 5 ) 

Taking 8 (n — 3)-dimensional combinatorial torus cube-tilings {w l, ^)x<j< 2 n-z with 1 < 
i < 8, one defines a n-dimensional rod ti/mg combinatorial cube packing 

(z\ w^') + [0, l] n for 1 < i < 8 and 1 < j < 2 n ~ 3 . 

Theorem 4.3. The probability of obtaining a rod tiling is 

Vf x q 8 n _ 3 

where q n is the probability of obtaining a n- dimensional cube-tiling andp\ 5 is a rational 
function of n. 

Proof. Up to equivalence, one can assume that in the random-cube packing process, 
one puts 

z 1 = (h 1 , w 1 ' 1 ) = fa, t 2l t 3 , . . . ) and z 2 = (h 2 , w 2 ' 1 ) = {h + 1, t 4 , h, . . . ). 

Then there are n(n — 1) possible choices for the next cube, 2(n — 1) of them are 
respecting the lamination. So, there are (n — 2)(n — 1) choices which do not respect 
the lamination and their probability is = Without loss of generality, we can 
assume that one has 

z 3 = (h 3 ^ 3 ' 1 ) = (t 6 ,t 2 + Lt 5 + !,•■• )■ 
In the next 5 stages we add cubes with n — 3 new parameters each. We have more than 
one type to consider under equivalence and we need to determine the total number 
of possibilities in order to compute the probabilities. 

For the cube z A + [0, l] n we should have three integers ii, i 2 , i 3 such that zf. = z?.+l 
(mod 2). Necessarily, the ij are all distinct, which gives n(n — l)(n — 2) possibilities. 
There are exactly 6 possibilities with ij < 3. One of them corresponds to the non- 
extensible cube packing of Figure [2] on the first 3 coordinates which the 5 others 
have a non-zero probability of being extended to the rod tiling. When computing 
later probabilities, we used the automorphism group of the existing configuration and 
gather the possibilities of extension into orbits. At the fourth stage, the 5 possibilities 
split into two orbits: 

(1) Of: (h\w^) for JG {1,2,3,4} with p\ = p 3 l n(n _ 3 {n _ 2) , 

(2) 0\: (h\w^) for i G {1,2,3,7} with p\ = p\ n(n J )(n _ 2) ; write = 3(n - 
3)(n — 4) + 3(n — 3) + 4 the number of possibilities of adding a cube to the 
packing ((h\w^) + [0,l} n ) ie{1X3J} . 
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When adding a fifth cube one finds the following cases up to equivalence: 

(1) 0\: (h\ w^ 1 ) for i G {1,2,3, 4, 5} with p\ = ^_2__ 

(2) 0\ 



(3) O 



(h\ w^ 1 ) for i 6 {1, 2, 3, 4, 7} with p\ = p\ 2(n J )(n _ 2) + pf ^ 
{h\ w*' 1 ) for i 6 {1,2,3, 7, 8} with p\ = p\^. 



When adding a sixth cube one finds the following cases up to equivalence: 



(1) 0\ 

(2) Q\ 

(3) 01 



(h l , w 

(h\w^) for i e {1,2,3,4,5,7} with V % = + vl^y 

(h\w^) for i e {1,2,3,4,7,8} with - - 5 1 1 - 5 3 



3(n- 


2) 


2 




3(n- 


2) 


1 




n(n- 


-2) 



J 3 3(n-2)- 

When adding a seventh cube one finds the following cases up to equivalence: 

(1) 0\: (ti, w*' 1 ) for i E {1,2,3, 4, 5, 6, 7} with p\ = p\ + ^_1_ 5 

(2) 0\: (h\ w*' 1 ) for i G{1,2,3, 4, 5, 7, 8} with p 2 7 = + p 6_^_ 

The combinatorial cube packing of eight cubes w i,1 ) + [0, l] n )i<i<8 is then obtained 
with probability p\ = p\ + pl^h^- 

Then we add cubes in dimension n — 4 following in fact the construction of Theorem 



4.2 The parameters £3, £4 and t 6 appear only two times in the cube packing for the 



rods, which contain 6 cubes in total. So, when one adds cubes we have 8(n — 3) choices 
respecting the cube packing, i.e. of the form z 9 = (h\w 1 ' 2 ) with w % - = w % - (mod 2) 
for some 1 < j ' < n — 3. We also have 3(n — 3)(n — 4) choices not respecting the rod 
tiling structure, i.e. of the form z 9 = (k 1 , w) with k l being one of h l for 1 < i < 3 with 
t 3 , t 4 or t 6 replaced by another parameter. But after adding a cube (h\ w h2 ) + [0, l] n 
with h % containing £3, £4 or t§ this phenomenon cannot occur. Below a type of 
probability p^ is a packing formed by the 8 vectors (h\w hl )i<i<s an d h — 8 vectors of 
the form (h l ,w h2 ) amongst which r of the parameters t 3 , £4 or do not occur. Note 
that there may be several non-equivalent cube packings with the same type but this 
is not important since they have the same numbers of possibilities. 
Adding 9 th cube one gets: 



(1) -T9 9 _ 8 2(n-3) 

V / J 3> Pi — Fl8(n-3)+3(n-3)(n-4)' 

(2) Tlpl^pl, ^ 



8(n-3)+3(n-3)(n-4) ■ 

Adding 10 i?i cube one gets: 

(1) n\pf=p\ — — 



7(n-3)+3(n-3)(n-4)> 

I'O^I TlO ^10 _ „9 6(n-3) . 9 3(n-3) 

l^J J 2 ' P2 - Pi 7( n _3) +3 ( n _3Vn-4l ^2 1 

(3) T/ , P3 = p\-, &=*L 



7(n-3)+3(n-3)(n-4) 1 ^2 7(n-3)+2(n-3)(n-4) ' 
) 4(n-3) 
!7(n-3)+2(n-3)(n-4) - 

Adding 11 th cube one gets: 

HI Til „!! _ y.10 6(n-3) , ip 2(n-3) 

V / J 2 ) Pi — Pi 6 (n~3)+3(n-3)(n-4) ^2 6(n-3)+2(ra-3)(ra-4) ' 

(n\ T n r, 11 — n 10 4(n-3) . 10 4(n-3) 

W J l > F2 —^2 6 ( n _3) +2 ( n _3)( n _4) "T- Pz 6 ( n _3) + ( n _3)( n _4); 

(Q\ T n n 11 — n 10 2(n-3) 

1°; J J ^3 — ^3 6(n-3)+2(n-3)(n-4)' 

Adding 12 th cube one gets: 

(i\ T 12 12 11 n^3 

V 1 ) L 2 ■> Pi — Pi 5(n-3)+2(ra-3)(ra-4)> 

T \2 12 _ 11 4(w-3) I „11 3(n-3) 

l z J J l iF2 — Pi 5(n-3)+2(n-3)(n-4) ^2 5(n-3)+2(n-3)(n-4) ' 



14 



MATHIEU DUTOUR SIKIRIC AND YOSHIAKI ITOH 



(Q\ T 12 n 12 — n 11 2(n-3) , n 5(n-3) 

V°) A iF3 ^2 5( n _3) +2 ( n -3)(n-4) ^ ^3 5(n-3)+2(n-3)(n-4) ' 

Adding 13 i?t cube one gets: 

m T 13 n 13 _ „12 4(n-3) , 12 2(n-3) 

l x J J l > Pi — Pi 4(n-3)+2(n-3)(n-4) ^2 4(n-3)+(n-3)(n-4) ' 

(n\ T 13 T) 13 - n 12 2(n-3) , i 2 4(n-3) 

l z J J ) —^2 4( n _3)+( n _ 3 )( n _4) "T- ^3 4 ( n _ 3 ) • 

Adding lA th cube one gets: 

(1) T 14 o 14 - » 13 

V / J l 'Pi Pi 3(n-3)+(n-3)(n-4)> 

fO^I T"14 „13 _ 13 2(n-3) . l 3 3(n-3) 

l z J J > Pi — Pi 3(„_3)+( n _3)( n _4) "^^2 3(n-3)' 

Adding 15 th cube one gets: 

fll T lb o 15 - o 14 2(n ~ 3) I- o 14 

l-U J J Pi — Pi 2(n-3)+(n-3)(n-4) ^^2 ■ 

After that if we add a cube z + [0, l] n , then necessarily z is of the form (h l ,w). 
So, we have 8 different (n — 3)-dimensional cube packing problems show up and the 
probability is pi 5 ^_ 3 - □ 
A combinatorial torus cube packing CV is called laminated if there exist a coordi- 
nate j and a parameter t such that for every cube z + [0, l] n of CP we have Zj = t 
(mod 1). 

Theorem 4.4. (%) The probability of obtaining a laminated combinatorial cube pack- 



mq is -. 



(ii) For any n>l, one has E(M£ (n)) < 2"(1 - + ^E(M^(n - 1)). 
fmj For any n > 3, ±E{Ml{n)) < 1 - fji 



Proof. Up to equivalence, we can assume that after the first two steps of the process, 
we have 

z 1 = (ti, . . . ,t n ) and z 2 = (t\ + l,t n +i, • • • ,*2n-i)- 
So, we consider lamination on the first coordinate. We then consider all possible 
cubes that can be added. Those cubes should have one coordinate differing by 1 with 
other vectors. This makes n(n — 1) possibilities. If a vector respects the lamination 
on the first coordinate then its first coordinate should be equal to t\ or t± + 1. This 
makes 2(n — 1) possibilities. So, the probability of having a family of cube respecting 
a lamination at the third step is -. But one sees easily that in all further steps, 
the choices breaking the lamination have a dimension strictly lower than the one 
respecting the lamination, so they do not occur and we get (i). 

By separating between laminated and non-laminated combinatorial torus cube 
packings, bounding the number of cubes of non-laminated combinatorial torus cube 
packings by 2 n one obtains 

E(ML(n)) < (1 - I) x 2" + l(E(M^(n - 1)) + E(Ml(n - 1))), 

which is (ii). (iii) follows by induction starting from |_E(M^(3)) = || (see Table 
§. □ 

5. Properties of non-extensible cube packings 

Theorem 5.1. If a combinatorial torus cube packing has at least 2 n — 3 cubes, then 
it is extensible. 
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Proof. Our proof closely follows [DIP 6] but is different from it. Take CV' a combina- 
torial torus cube packing with 2 n — a cubes, a < 3. Take N such that Nb(CV, N) > 
and CV a discrete cube packing with 4>{CV) = CV' . If CV is extensible then CV' is 
extensible as well. 

We select 8 £ R and denote by Ij the interval [S + i,5+ for < j < 3. Denote 
by rij t k the number of cubes, whose k-th coordinate modulo 2 belong to Ij. 

All cubes of CV, whose k-th coordinate belongs to Ij, Ij + \ form after removal of 
their k-th. coordinate a cube packing of dimension n — 1, which we denote by CVj^- 
We write rij^ + = 2™ _1 — dj k and obtain the equations 



We can then write the vector dk = (c?o,fc, di,k, <^2,fc, cfa.fc) in the following way: 



4 = ci(l,l,0,0) + c 2 (0,l,l,0) + c 3 (0,0,l,l) + c 4 (l,0,0,l) with J2 c i= a 



and Q G Z + . This implies = Cj + Cj+i < ^2 Cj = a. This means that the (n — 1)- 
dimensional cube packing CVj^ has at least 2 n_1 — 3 cubes, so by an induction 
argument, we conclude that CVj t k is extensible. 

Suppose now that the k-th coordinate of the cubes in CV have values < 5i < 62 < 
■ ■ ■ < 5m < 2. So, the set of points in the complement of CV, whose k-th coordinate 
belongs to the interval [Si, S i+ i[ with 5m+i = 5i + 2 can be filled by translates of the 
parallelepiped Paral k (a) = [0, l[ fc_1 x [0, a{x [0, l[ n ~ k . 

Note that as 5 varies, the vector dk varies as well. Suppose that for some i, we 
have the k-th layer [Si, <5 i+1 [ being full and [<5j_i, Si[ containing x translates with x < 3 
of the parallelepiped Paralk(Si + \ — Si). Then if one selects another coordinate k', 
all parallelepipeds Paralk'(S' i/+1 — S[,) filling the hole delimited by the parallelepiped 
Paralk{Si + \ — Si) will have the same position in the k-th coordinate. This means that 
they will form x cubes and that the cube packing is extensible. This argument solves 
the case a = 1, because up to symmetry d k = (0, 1, 1, 0). 

If a = 2, then the case of vector of coordinate dk being equal to symmetry to 
(0, 2, 2, 0) or (0, 1, 2, 1) is also solved because we have seen that a full layer implies 
that we can fill the hole. We have the remaining case (1, 1, 1, 1). If the hole of this 
cube packing cannot be filled, then we have a structure of this form: 



Selecting another coordinate k', we get that the two parallelepipeds z + Paralk{Si + \ — 
Si) and z' + Paralk{Si' + i — S^) have zi = z[ for / 7^ k, k! . This is impossible if n > 4. 
So, if dk = (1, 1, 1, 1) for some k and S, then the hole can be filled. 



3 





j=0 



4 
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If a = 3, and d k , up to symmetry, is equal to (0, 3, 3, 0) or (0, 2, 3, 1) then we have 
a full layer and so we can fill the hole. If the vector d k = (2, 1, 1, 2) occurs, then by 
the same argument as for (1, 1, 1, 1) we can fill the hole. □ 

Proposition 5.2. (i) Non extensible combinatorial torus cube packings of dimension 
n have at least n + 1 cubes. 

(ii) If CV is a combinatorial torus cube tiling, then in a coordinate j a parameter 
t occur the same number of times as t and t + 1. 

Proof, (i) Suppose that a combinatorial torus cube packing CV has m < n cubes 
(z l + [0, l] n )i<i< m . By fixing Zi = z\ + 1 for % = 1,2, ... ,m we get that the cube 
z + [0, l] n does not overlap with CV. 

(ii) Without loss of generality, we can assume that a given parameter t occurs only 
in one coordinate k as t and t + 1 . The cubes occurring in the layer [t,t+l], [t + l,t + 2] 
on j-th coordinate are the ones with Xj = t, t + 1; we denote by V t and V t+ i their 
volume. Now if we interchange t and t + 1 we still obtain a tiling, so V t < V t+ i and 
V t+ i < V t . So, V t = V t+ i and the number of cubes with Xj = t is equal to the number 
of cubes with Xj — t + 1. □ 

Take a combinatorial torus cube packing CV obtained with strictly positive prob- 
ability. Let us choose a path p to obtain CV. Denote by N k , p (CV) the number of 
cubes obtained with k new parameters along the path p. 

Proposition 5.3. Let CV be a non- extensible combinatorial torus cube packing, p a 
path with p(CV,p, oo) > 0. 

(i) N nyP {CV) = 1 and N n _ hp (CV) = 1. 

(ii) N n -2, P (CV) < 2 and N n _2, P (CV) = 2 if and only if CV is laminated, 
(in) One has N kp (CV) > 1 for < k < n. 

M N(CV) = Y2 =0 kN Kp {CV) > 

(v) IfN(CV) = T -^f^ then N Kp (CV) = 1 for k > 1. 

Proof. The first cube z 1 + [0, l] n has n new parameter, but the second cube z 2 + [0, 1]™ 
should not overlap with the first one so it has n — 1 parameters and N n ^ p {CV) = 1. 
Without loss of generality, we can assume that z 1 = (t 1: ... : t n ) and z 2 = (ti + 
1, t n+ i, . . . , i2n-i)- When adding the third cube z 3 + [0, l] n , we have to set up 2 
coordinates depending on the parameters ti, i < In — 1 thus N n _i p (CV) = 1. 

If z\ — t\ or t\ + 1 then we have a laminated cube packing, we can add a cube with 
n — 2 parameters and N n _2, P (CV) = 2. Otherwise, we do not have a laminated cube 
packing, three coordinates of z 4 need to be expressed in terms of preceding cubes and 
thus N n _ 2yP {CV) = 1. 

(iii) The proof is by induction; suppose one has put m' = Y^i=k Ni, P (CV) cubes. 
Then the cube z m ' + [0, l] n has k new parameters t\, . . . , t' k in coordinates i 1: . . . , i k . 
The cube C = z + [0, l] n with z ix — t[ + 1 and z { = zf for i ^ . . . , i k } has — 1 
free coordinates {i 2 , ■ ■ ■ ,i k } and thus k — 1 new parameters. So, N k _i jP (CV) > 1. 

(iv) and (v) are elementary. □ 

Conjecture 5.4. Let CV be a combinatorial torus cube packing and p a path with 
p(CV,p, oo) > 0. 

(i) For all k > 1 one has ^' fc=0 N n _ k>p (CV) < 2 l . 
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(ii) N{CV) < 2 n - 1; if N(CV) = 2 n - 1, then CP is obtained via a lamination 
construction. 



A "perfect matching of a graph G is a set M. of edges such that every vertex of G 
belongs to exactly one edge of Ai. A 1- factorization of a graph G is a set of perfect 
matchings, which partitions the edge set of G. The graph K 4 has one 1-factorization; 
the graph K 6 has, up to isomorphism, exactly one 1-factorization with symmetry 
group Sym(5). 



Proposition 5.5. Let CV be a non- extensible combinatorial torus cube packing. 

(i) If n is even then CV has at least n + 2 cubes. 

(ii) Ifn is odd andCV has n+1 cubes then N(CV) = n( - n ^ . Fix a coordinate j and 
a parameter t occurring in at least one cube. Then the number of cubes containing t, 
respectively t + 1 in coordinate j is exactly 1. 

(Hi) Ifn is odd then isomorphism classes of non- extensible combinatorial torus cube 
packings with n + 1 cubes are in one to one correspondence with isomorphism classes 
of 1-factorizations of K n+1 . 

(iv) If n is odd then the non- extensible combinatorial torus cube packings with n + 1 
cubes are obtained with strictly positive probability and f^(n) = f> Qoo (n) — n+1. 



Proof. We take a non-extensible cube packing CV with n + 1 cubes. Suppose that for 
a coordinate j we have two cubes z l + [0, l] n and z l + [0, l] n with Zj = z* = t. If a 
vector z has Zj = t + 1, then z + [0, l] n does not overlap with z 4 + [0, l] n and z l +[0, l] n . 
There are n — 1 remaining cubes to which z + [0, l] n should not overlap but we have 
n — 1 remaining coordinates so it is possible to choose the coordinates of z so that 
z+ [0, l] n does not overlap with CV. This is impossible, therefore parameters appear 
always at most 1 time as t and at most one time as t + 1 in a given coordinate. 

By Lemma 4.1 every parameter t appear also as t + 1. So, every parameter t appears 
one time as t and one time as t + 1. This implies that we have an even number of 
cubes and so (i). Every coordinate has ! ^ parameters, which gives n( - n 2 +1 - > parameters 
and so (ii). 

(iii) Assertion (ii) implies that any two cubes Ci and CV of CV have exactly one 
coordinate on which they differ by 1. So, every coordinate correspond to a perfect 
matching and the set of n coordinates to the 1-factorization. 

(iv) Since parameters t appear only one time as t and t+1, the dimension of choices 
after k cubes are put is n — k and one sees that such a cube packing is obtained with 
strictly positive probability. The existence of 1-factorization of K 2p (see, for example, 
[A106, Ha78j) gives f^in) < /> 0o o( n ) — n + 1- Combined with Theorem 5.2 i, we 
have the result. □ 



Conjecture 5.6. Ifn is even then there exist non- extensible combinatorial torus cube 
packings with n + 2 cubes and w( " 2 +1 ^ parameters. 
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In dimension 4 there is a unique cube packing (obtained with probability ^) 
satisfying this conjecture: 



/ 


tl 


*2 


*3 


u \ 




h 


*6 


h 


U + l 




h + 1 


h 


t 7 + l 


tg 




£ 5 + i 


h + i 


t 3 + l 


tw 




ti + 1 


h + i 


t 7 


ho + 1 


V 


h 


t 2 + 1 


t 7 + l 


t 9 + l J 



Proposition 5.7. (i) There are 9 isomorphism types of non- extensible combinatorial 
torus cube packings in dimension 6 with 8 cubes and at least 21 parameters (see Figure 
; they are not obtained with strictly positive probability. 
W8 = /£(6)<# 0)DO (6). 

Proof, (ii) follows immediately from (i). The enumeration problem in (i) is solved 
in the following way: instead of adding cube after cube like in the random cube 
packing process, we add coordinate after coordinate in all possible ways and reduce by 
isomorphism. The computation returns the listed combinatorial torus cube packings. 
Given a combinatorial cube packing CP in order to prove that p{CV, 00) = 0, we 
consider all (8!) possible paths p and see that for all of them p(CV,p, 00) =0. □ 

Proposition 5.8. Ifn — 3, 5, 7, 9, then there exist a combinatorial torus cube tiling 
obtained with strictly positive probability and n ^" 2 +1 - > parameters. 

Proof. If n is odd consider the matrix H n = rriij with all elements satisfying rrii+k,i = 
rrii_k,i + 1 for 1 < k < the addition being modulo n. The matrix for n = 5 is 



h 


t 7 + l 


tl 3 + 1 


tu 


tio 




*2 


*8 + l 


t u + l 


tl5 


tu 


t7 


h 


tg + 1 


tl 5 + 1 


iii + l 


tl2 


ts 


u 


h + i 


*6 + l 


£l2 + l 


tl3 


h 


h 



Then form the combinatorial cube packing with the cubes (z l + [0, l] n )i<i<„ and z % 
being the z-th row of H n . It is easy to see that the number of parameters of cubes, 
which we can add after z l is n— i. So, those first n cubes are attained with the minimal 
number n(n 2 +1 ^ of parameters and with strictly positive probability. If a cube z+[0, l] n 
is non-overlapping with z l + [0, l] n for i < n then there exist a(i) G {1, . . . ,n} such 
that z a (i) = z*,^ + 1. If i ^ i' then a(i) 7^ o~(i'), which proves that a G Sym(n). 
We also add the n cubes corresponding to the matrix H n + Id n . So, there are n! 
possibilities for adding new cubes and we need to prove that we can select 2 n — In 
non-overlapping cubes amongst them. 

The symmetry group of the n cubes (z l + [0, l] n )i<j<„ is the dihedral group D 2n 
with In elements. It acts on Sym(n) by conjugation and so we simply need to list 
the relevant set of inequivalent permutations in order to describe the corresponding 
cube packings. See Table [3] for the found permutation for n — 3, 5, 7, 9. □ 

The cube packing of above theorem was obtained for n = 5 by random method, 
i.e., adding cube whenever possible by choosing at random. Then the packings for 
n = 7 and 9 were built using the matrix H n and consideration of all possibilities 
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( ti ts tg tu + 1 tir + 1 tig ^ 
tl + 1 <6 tlO tl3 + 1 tl6 + 1 tig 

t2 ts + 1 ill tl3 tl8 tig 

ti + 1 <7 tg + 1 iis ti6 i2i 

*3 <6 + 1 tl2 tl4 tlS + 1 <21 + 1 

t 3 + l i8 *10 + 1*16 + 1 *17 *20 + l 
*4 *7 + 1 tl2 + 1 tl3 + 1 tl7 + 1 tig + 1 
\t 4 + 1 t 8 + 1 ^11 + 1 tl4 + 1 tie + 1 tig + 1/ 

21 parameters, \Aut\ = 4 



/ ti is tg 

11 + 1 i6 iio 
ti is + 1 iio + 1 

12 + 1 ie + 1 tg + 1 



tlS 

ii4 
iis 
iie 

tis + 1 iis + 1 i 22 + 1 



ii7 


t21 \ 


tl8 


til 


ii9 


i22 


i20 


i22 



tn 


tii \ 


tin 


tii 


tig 


tii 


i20 


tii 



tz t-j tn 

is + 1 is ii2 ii4 + 1 ii? + 1 tn + 1 
U t 7 + 1 ii2 + 1 tie + 1 iig + 1 tn + 1 
\i 4 + 1 is + 1 in + 1 tie + 1 i2o + 1 hi + l/ 
22 parameters, \Aut\ = 64 

/ ti is tg ti3 
ii + 1 i6 iio fi4 
*2 is + 1 iio + 1 iis 
ti + 1 i 6 + 1 t 9 + 1 tie 

ts i? tn tis + 1 1 20 + 1 i 22 + 1 
ts + 1 is ti2 tie + 1 ii? + 1 i22 + 1 
U t 7 + l tu + 1 ti4 + 1 tig + 1 t 2 i + 1 
\t 4 + 1 is + 1 in + 1 iis + 1 tis + 1 t2i + lj 
22 parameters, \Aut\ = 16 

f ti ts tg ti3 tl7 t21 \ 

ti + 1 ie tio ti4 tis i22 

t2 is + 1 in iis tl9 i22 + 1 
t2 + 1 ie + 1 tg + 1 ti6 i20 i2i 

t3 t7 tlO + 1 tl3 + 1 t20 + 1 t22 

t 3 + 1 t 8 tl2 tl4 + 1 tig + 1 t21 + 1 
U t 7 + 1 tl2 + 1 tl6 + 1 tl8 + 1 t21 + 1 

\t 4 + 1 is + 1 in + 1 tie + 1 ti7 + 1 i22 + 1/ 
22 parameters, \Aut\ = 32 

/ ti ts tg 

tl + 1 t6 tlO 

ti is + 1 in 
t2 + 1 ie + 1 tg + 1 



tl7 


til \ 


iis 


til 


tl9 


til 


tio 


til 



( ti tg tg ti3 

ti + 1 te tio ti4 
ti is + 1 iio + 1 iis 
ti + 1 1 6 + 1 tg + 1 tie 

t 3 t 7 tn tlS + ltl8 + lt22 + l 

ts + 1 is tl2 tl4 + 1 tl7 + 1 t 2 2 + 1 
t4 t7 + 1 tl2 + 1 tlS + 1 t20 + 1 t 2 l + 1 

\t 4 + 1 is + 1 in + 1 tie + 1 tig + 1 t 2 i + 1/ 
22 parameters, \Aut\ — 64 



/ ti ts i9 

tl + 1 t(j ti() 

ti is + 1 tio + 1 
t2 + 1 te + 1 tg + 1 

t3 t7 til 

t3 + 1 ts tl2 



ii3 ii7 i2i \ 
tl4 iis t22 
iis ii9 i22 

il6 t20 t21 

tis + 1 1 20 + 1 hi + 1 

tl4 + 1 tl9 + 1 *21 + 1 



U t 7 + 1 tl2 + 1 tl6 + 1 tl7 + 1 t22 + 1 
\t 4 + 1 t 8 + 1 tn + 1 tlS + 1 tl 8 + 1 t 2 l + 1/ 

22 parameters, \Aut\ = 16 



\ 



/ tl ts t9 tl3 tl7 t21 

ti + 1 te iio ii4 tis i2i 
i2 is + 1 tn iis iis + 1 i22 
i2 + 1 te + 1 tg + 1 tie ii9 i22 

t3 t7 tlO + 1 tl3 + 1 t20 t22 + 1 

ts + 1 t8 tl2 tl4 + 1 tl7 + 1 t22 + 1 

t 4 t 7 + 1 tl2 + 1 tlS + 1 tl9 + 1 t 2 l + 1 

\t 4 + 1 1 8 + 1 tn + 1 tie + 1 120 + 1 t 2 i + 1/ 
22 parameters, \Aut\ — 16 



\ 



/ ti ts t9 tl3 tl7 t21 
ti + 1 te iio ii4 tis i2i 
i2 is + 1 tn iis iis + 1 i22 

t2 + 1 te + 1 tg + 1 ti6 tl9 t22 

t3 t7 tio + 1 tl3 + 1 t20 t22 + 1 

t 3 + l t8 tl2 tl 5 + ltl9 + lt21 + l 

t4 *7 + 1 tl2 + 1 tl4 + 1 tl7 + 1 t22 + 1 

\t 4 + 1 1 8 + 1 tn + 1 tie + 1 t2o + 1 t 2 i + 1/ 
22 parameters, \Aut\ — 8 



tis 

tl4 
tlS 
tl6 



tl7 t21 
tlS t22 

tl9 t22 + 1 
t20 t21 



t3 t7 tlO + 1 tl3 + 1 t20 + 1 t22 

i3 + 1 tg tl2 tl5 + 1 tlS + 1 t21 + 1 

U t 7 + 1 tl2 + 1 tl4 + 1 tig + 1 t21 + 1 
\t 4 + 1 t 8 + 1 tn + 1 tie + 1 tl7 + 1 t22 + 1/ 

22 parameters, \Aut\ = 16 



Figure 3. The non-extensible 6- dimensional combinatorial cube pack- 
ings with 8 cubes and at least 21 parameters 



invariant under the dihedral group L>2n by computer. But for n — 11 this method 
does not work. It would be interesting to know in which dimensions n combinatorial 
torus cube tilings with ra( -" 2 +1 - ) parameters do exist. 
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in — \ 


(1,2,3) 


n = 5 


(1,2,3,4,5) 


(1,2)(3,5,4) 


(1,4,5,3,2) 


n = 7 


(1,2,3,4,5,6,7) 
(1,7)(2,5,6)(3,4) 
(1,3,7)(2,6)(4,5) 


(1,7)(2,5,4,3,6) 
(1,6,2,3,7)(4,5) 
(1,3,7)(2,5,4,6) 


(1,6,2,5,4,3,7) 
(1,7)(2,3,4,5,6) 
(1,5,4,3,2,6,7) 


n = 9 


(1,2,3,4,5,6,7,8,9) 
(1,5,9)(2,8)(3,6,7,4) 
(1,5,4,3,9)(2,8)(6,7) 

(1,6,5,4,7,8,2,3,9) 
(1,9J(2,3,4, r,6,5,8J 
(1,6,2,3,9)(4, 7,8,5) 
(1,7,6,4,9)(2,5)(3,8) 
(1,4,9)(2,5)(3,7,6,8) 
(1,4,5,2,8,9)(3,7,6) 

(1,7,6,5,4,3,2,8,9) 


(1,6,7,4,3,5,9)(2,8) 

(1,9)(2,5,4,3,6,7,8) 

(1,9)(2,5,6,7,8)(3,4) 

(1,9)(2,3,6,7,8)(4,5) 

(1,9)(2,3,4, 5,6, 

(1,5,4,7,8,6,2,3,9) 

(1,7,6,9)(2,8,3,4,5) 

(1,4,8,3,7,6,9)(2,5) 

(1,4,9)(2,5)(3,8,7,6) 

(1,9)(2,5,8)(3,4)(6,7) 


(1,5,6,7,4,3,9)(2,8) 

(1,6,7,8,2,5,4,3,9) 

(1,9)(2,3,6,5,4, 7,8) 

(1,6,7,8,2,3,9)(4,5) 

(1, 9J(2, 3, 6J(4, l,6,b) 

(1,8,3,7,6,4,9)(2,5) 

(1,4,5,2,8,3,7,6,9) 

(1,7,6,3,4,5,2,8,9) 

(1,3,7,6,9)(2,8,4,5) 



Table 3. List of permutation describing combinatorial torus cube 
tilings with parameters in dimension 3, 5, 7, 9 
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